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Stochastic webs were discovered, first by Arnold for multi-dimensional Hamiltonian systems, and later by 
Chernikov et al. for the low-dimensional case. Generated by weak perturbations, they consist of thread-like 
regions of chaotic dynamics in phase space. Their importance is that, in principle, they enable transport from 
small energies to high energies. In this introductory review, we concentrate on low-dimensional stochastic webs 
and on their applications to quantum transport in semiconductor superlattices subject to electric and magnetic 
fields. We also describe a recently-suggested modification of the stochastic web to enhance chaotic transport 
through it and we discuss its possible applications to superlattices. 
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1 Introduction 

Stochastic webs exist in the phase spaces of Hamiltonian systems, that is, in the space formed 
by the coordinates and momenta of a dynamical system evolving in the absence of dissipation. 
They consist of a network of very thin thread-like regions within which the dynamics is chaotic, 
whereas the dynamics remains regular everywhere else. Although the concept seems abstract 
and mathematical at first sight, stochastic webs are now known to arise in a number of practical 
contexts, including for example plasma physics [H, ultra-cold atoms in optical lattices (3,0,0] 
and electrons in semiconductor superlattices (SLs) [1, 0, 0, 0, 0, 0, U]; they have also been 
considered in connection with celestial mechanics [121 ]. The importance of the chaotic threads 
is that they can transport matter and energy effectively over long distances 0,113. In this 
brief and rather informal review we aim to introduce the general reader to stochastic webs, 
explaining what they are and discussing their recent developments and applications, taking 
electron transport in semiconductor SLs as our example. 

We start (Section 11.11) from the definition of a Hamiltonian system, its dimensionality and 
integrability. Then (Section [L2]) we consider the effect of perturbations of the integrable system, 
which brings us to the concept of a chaotic (stochastic) layer, in particular related to resonances. 
The latter allows us to explain in the beginning of Section [2] the onset of the Arnold stochastic 
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web in multi-dimensional systems. The main purpose of this section is to discuss the more 
sophisticated nature of low-dimensional stochastic webs which are, however, still related to the 
concept of resonance. In Section [3J we explain the limitations of the transport though the low- 
dimensional web and suggest a subtle way of overcoming these limitations. Finally, in Section [4j 
we discuss a rather unexpected application of the stochastic web concept to quantum electron 
transport in nanometer-scale semiconductor SLs in the presence of electric and magnetic fields. 
Section [5] draws conclusions. 



1.1 Hamiltonian systems 

Hamiltonian systems play an important role in physics, chemistry, biology and engineering, and 
form a fundamental class of dynamical systems 0, [TtJ ■ They are defined by the following 
dynamical equations: 



dp^ _ _dH_ dqi L _dH_ 
dt dqi ' dt dpi ' 

If the Hamiltonian H does not depend on time t, while depending only on the momenta p 
{ p ..... /) \ } and coordinates q = (qi, . . . ,qw), then it is called ^-dimensional. If it also depends 
on time t, then it has the dimension N + |. A remarkable property of any Hamiltonian system 
is the equality of its full and partial derivatives with respect to time: 



dH dH 

In particular, for time-independent Hamiltonians, H(p, q) is conserved along the trajectory. 
In general, the equations of motion (1) may not be integrable in quadrature^] (lEI . 16 . 17]. 



whence the importance of integrable systems, i.e. those time-independent Hamiltonian systems 
for which a transformation {p, q} <-> {/, 9} exists such that 

H(p,q)=H(l). (3) 

Ii are called actions while 6i are called angles. It follows from (3) that I is conserved: 

^ = -^ = (4) 
dt dOi ' 1 ' 

while the angles 0, change with constant speeds (for a given I), 

d6i _ dH{T) 



dt dL 



(5) 



which are called frequencies. 

Note that angles 9{ are cyclic variables i.e. p and q are periodic functions of 9i with a period 
2ir for any 9{ |l|, [l5|, 0, [l3]. Thus, Eqs. (4) and (5) correspond to periodic or quasi-periodic 
motion. The simplest and most often used example of an integrable system is a one-dimensional 
one, which will be described in more detail below. 



1 When the solution of a differential equation expressible in terms of a formula involving integrations, it is said to be solvable 
by quadrature. 
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Figure 1. Schematic diagram to show a weak distortion of the majority of trajectories by a weak time-periodic perturbation: 
the blue line shows the trajectory of the unperturbed system, while the red line shows the stroboscopic Poincare section of 
the trajectory of the perturbed one. 




-71 71 - x 

Coordinate q (mod 2ti) Coordinate q (mod 2u) 



Figure 2. (a). The separatrix of the pendulum H = Ho = p /2 — cos(q): the separatrix corresponds to H = H s = 1. (b). 
The chaotic layer (replacing the separatrix) in Poincare section of the ac-driven pendulum: H = Hq — O.Olqcos(t). 



1.2 Perturbed Hamiltonian systems 

It is natural to pose the question: what is the effect of a weak perturbation on an integrable 
system? For the majority of cases, the answer is given by the Kolmogorov-Arnold-Moser (KAM) 



theory 15l |: most of the trajectories are just weakly distorted by a weak perturbation while 
remaining regular. Let us illustrate this by an example of a one-dimensional system weakly 
perturbed time-periodically. In this case, it is convenient to present the trajectory in the stro- 
boscopic Poincare section [l|, 15, 16, [13], i.e. presenting states of the system {p(t),q(t)} only 
at the discrete sequence of instants t = t n = t$ + nT where to is some initial instant, T is the 
perturbation period and n = 0, 1, 2, .... If the trajectory is "just weakly distorted while remain- 
ing regular", then, in particular, the unperturbed trajectory and the Poincare section of the 
perturbed one have the same dimension of 1 (i.e. they are just lines), and the same topology, 
while just slightly deviating from each other (Fig. [IJ). 

There are, however, two kinds of situation for which KAM-theory is not valid. The first of 
these relates to the separatrices of the unperturbed systems. Let an integrable system possess a 
separatrix i.e. the line (or surface, or hyper-surface in the general multi-dimensional case) that 
separates trajectories of a different topology in the phase spac^J e.g. in the example shown 
in Fig. Eta), the separatrix separates closed trajectories (corresponding to oscillations inside 
the separatrix loops) from open trajectories (corresponding to the running coordinate below or 
above the separatrix loops). If the system is perturbed time-periodicalljH, then the separatrix is 
replaced by a chaotic trajectory. In Poincare section, the chaotic trajectory lies within a chaotic 
layer (Fig. [2(b)): the latter has a complicated (fractal) structure but its outer boundaries are 
well defined and the region delineated by these boundaries has the dimension 2, unlike the 



1 More rigorously, the separatrices may be defined as follows I18H . Let the integrable system possess a saddle i.e. a hyperbolic 
point in the one-dimensional case (i.e. an unstable stationary point with an exponential dynamics of trajectories approaching 
it), or a hyperbolic invariant torus in higher-dimensional cases. The stable (incoming) and unstable (outgoing) manifolds 
are called separatrices. 

2 In multi-dimensional cases, a time-independent perturbation also may give rise to the invalidity of the KAM-theory near 
the separatrix. 
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E (energy) 

Figure 3. A schematic diagram showing the dependence of eigenfrequency lu on the energy E of an eigenoscillation, and 
the meaning of the resonance energy E r . 

dimension 1 of regular trajectories. Thus, even the appearance of the Poincare section allows us 
to distinguish immediately between regular and chaotic trajectories, unless of course the width of 
the chaotic layer is less than the accuracy provided by the numerical integration of equations of 
motion. The theoretical prediction of the width in energy of the chaotic layer has a long and rich 
history. Its description on a physics level of rigour may be found in the book by Zaslavsky (l7| . 
Studies that are more mathematically rigorous have recently been reviewed [19(. The maximum 
width of the layer and other significant features (hi gh p eaks) of the width as function of the 
perturbation frequency have recently been described [20l. I21I] . 

Another characteristic situation where the KAM-theory is invalid relates to resonances, i.e. to 
areas of the phase space where at least one of the following conditions holds 



noJi{I r ) = mu)j(I r ), (6) 
n,m = ±1, ±2, ±3, . . . , 

3 = 1,2,3, ... ,iV, i^j, 
AT = 2,3,4,..., 



or 



noJi(I r ) =luf, (7) 
n = ±1,±2,±3, . . . , 

i = 1,2,3,... ,iV- i, 

N= 3 -, 5 -, 7 -,..., 
2'2'2' ' 

where uif is the frequency of the corresponding time-periodic perturbation while I is the number 
of the Fourier harmonic that may exist for the time-periodic perturbation (e.g. for a monochro- 
matic perturbation, only I = 1 is relevant). 

The rigorous treatment of motion in the resonance range is rather complicated, being related 
to the Poincare-Birkhoff theorem and homoclinic and heteroclinic trajectories and tangencies 
12) 23[. We do not consider it here. Rather, we give a brief interpretation of the resonance-related 



chaos in physical terms^- For the sake of clarity, consider an ac-driven ID Hamiltonian system 
whose frequency of eigenoscillation u increases monotonically with the energy of eigenoscillation 
E = Ho(p,q), while the perturbation frequency LVf exceeds the minimum of oj{E), as shown in 
Fig. El 



This was given for the first time by Chirikov [24ll (a clear presentation of the issues in question can be found e.g. in |17|). 
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H = H (p,q) - hqcos(u f t), (8) 
h < 1, 

ujf > luq = mm{uj(E)}. 

Then there necessarily exists an energy E r such that the resonance condition (7) with n = I = 1 
holds true: 

u{E r )=oj f . (9) 

Consider motion for energies close to E r . Let us transform from variables {p, q} to action-angle 
variables {1,9}, so that the Hamiltonian becomes 

H(p, q) = H(I, 6)= I' dlu(l) -h^ Qn(I) cos(n9) cos(uJ f t) 

J Imi„ n 

^ ~ ~ l 

= / dlco(I) - -hqi cos(9 — ujjt) + ... 
It £ 

= H (I,9 = e-u; f t) + V f (I,e,t), (10) 
I{p,q) = 1(E) = i- j> p{q,E) dq, E = H (p,q), 

1 f 27T 

q n = q n (I) = ^ J d9qcos(n9). 

Here, the dots "..."denote terms that vary with time much faster than the preceding term: 
they are denoted in the next equality as Vf. Thus, allowing for the resonance condition (9), we 
may introduce the slow angle 9 = 9 — to ft and present the original Hamiltonian as a sum of an 
"autonomous" part Hq(I,9) and the time-dependent (fast-oscillating) part Vf(I,9,t). It is easy 
to check by direct substitution into the corresponding Hamiltonian equations of motion that the 
dynamics of the variables I and 9 is governed by the Hamiltonian 

H = H-io f I. (11) 

Taking into account that the perturbation is small (h ^ 1) and, therefore, that the variation of 
/ around the resonance value I r = I(E r ) is also small, we may approximate the function uj(I) 
near the resonance value as 



U(I) W UJf +u' r (I- I T ), 



, du 



I=I r 



(12) 



From (10)-(12), we ultimately obtain the approximate auxiliary Hamiltonian governing the dy- 



November 11, 2009 16:31 Contemporary Physics CPstochasticwebs'rcsub 



6 Soskin, McClintock, Fromhold, Khovanov & Mannella 

namics of {I, 9}: 

8(1, 0, t) = l -J r (l - I r ) 2 - l -h Ql cos(0) + Vf(I, 9, t) 

= H (I = I-I r ,9) + V f . (13) 

It represents the sum of a pendulum- lika3 autonomous Hamiltonian Hq(I, 9) and a time-periodic 

(fast-oscillating) part Vf that plays the role of a perturbation. The pendulum-like part Hq 
possesses a separatrix (cf. Fig. [2(a)) while the perturbation- like part Vf tends to destroy the 
separatrix, replacing it with an exponentially narrow chaotic layer. 

Thus we have shown that a resonance is necessarily associated with a narrow chaotic layer. 



2 Stochastic webs 



In the example considered above, the chaotic layer associated with the resonance provides only 
a narrow (oc yh) variation of energy (or, equivalently, of action). Thus, there is no significant 
transport in energy. Let us pose a question: could there be situations when a perturbation provides 
for chaotic transport through a large range of energies? We now describe the three stages of 
conceptual evolution that led to a positive answer to this question. 



2.1 Multi- dimensional web 



First, Arnold showed in 1964 25f] through rather simple topological arguments (also presented 
clearly in [13]) that, if the system is multi-dimensional (namely, if N > 5/2), and if the so called 
non-degeneracy condition det(d 2 H/dlidlj) ^ is fulfilled (in other words, if the system is 
sufficiently nonlinear), then resonances necessarily intersect with each other, forming an infinite 
web in the phase space along which an exponentially slow chaotic diffusion occurs. 



2.2 Low-dimensional webs 

Secondly, Chernikov et al. published an important series of papers in the late 1980s. We shall 
review just three of the more important of them, concentrating on the model of a harmonic 
oscillator subject to a plane wave, which will be relevant to our discussion of semiconductor SLs 
below. A good review of a broad range of the early work on low-dimensional stochastic webs 
may be found in [J]; more recent work is reviewed in [26| (see also [l7|). 

The main idea of Chernikov et al. is that a stochastic web may arise even in low-dimensional 
systems (N = |;2) provided that the non-degeneracy condition is lifted, in other words, in this 
case, if 

duj 

while the perturbation in the equation of motion is resonant and coordinate-dependent. 
2.2.1 Cobweb 

We now consider the best known example of a low-dimensional stochastic web, the skeleton of 
which in p — q Poincare section has a form resembling that of a cobweb (Figs. 13(a), E( D ))- We 



(14) 



1 I plays the role of a generalized velocity while 9 plays the role of the generalized coordinate. Note that the generalized 
potential contains the small multiplier h, so that the maximal absolute value of the "velocity"/ is proportional to yh and, 
therefore, is small too. 



November 11, 2009 16:31 Contemporary Physics CPstochasticwebs'rcsub 

Stochastic webs and quantum transport 



P 
0.6 



0.0 



-0.6 





(a) 
























-0.1 



0.0 



0.1 x/2it 



Figure 4. Stochastic web for the system (15), as obtained numerically with cjo = 1, v = 4, k = 15, euj^/k = 0.1 [27ll . The 
integration time is 2x 10 4 periods of oscillation 2n/uo- (a). Poincare section, (b). The corresponding probability distribution. 

suppose that a harmonic oscillator is perturbecd by a resonant plane wave 27]: 



q + aigf? = e-77 sin(feg — ut), 



(15) 



v = nujQ, 



n 



1,2,3, 



This particular model has a number of applications, especially to plasma physics [1[ and to 
semiconductor SLs, as shown below in Sec. HI In order to understand the origin of the stochastic 
web shown in Fig. HI we 



(i) transform to polar coordinates {p, 6} or, equivalently, to action-angle variables {1,9}: 

q = psm(9), p = q = u>opcos(9), 

P 



(16) 



(ii) make use of the formula [2 



cos(xsin(0) — y) = J m (x) cos(m9 — y) 



(17) 



m=0 



where J m (x) is a Bessel function of the mth ordeiOJ 



Using (16) and (17), it is not difficult to show that the Hamiltonian of a harmonic oscillator 
perturbed by a plane wave can be represented in action-angle variables as 



H(I, 9, t) = lo I + e ^Yl J m(kp(I)) cos{m9 - ut). 



(18) 



Note that, due to the resonance condition v = uloq in (15), the term in the sum in (18) cor- 
responding to m = n is nearly constant compared to other terms in the sum. So, similarly to 



2 For parameters, we use the same notation as Zaslavsky p] and Chernikov et al. 1271 while the coordinate is denoted as q 
(instead of x in Q]. 2V : cf. Fig. 4(a)) in order to match the notation in other sections and in some figures from other works 
reproduced below. 

Note that J m (x) is an oscillatory function of x with gradually decreasing amplitude as x increases. At X ~ 1, the period 
of oscillation is ~ 2tt while the amplitude is ~ 1. For large x, the Bessel function asymptotically approaches the function 
yiTT^ra) cos(z - (2m + 1)tt/ 4 )- 
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Figure 5. Left panel: schematic representation of the grid-like separatrix of the Hamiltonian H s , as defined in Eqs. (20); 
saddles are shown by dots. Right panel: schematic representation of the same separatrix in Poincare section p — q. 

the case of nonlinear resonance considered in Section Tl.21 above, it is this term that provides the 
major contribution to the dynamics. The other terms in the sum play the role of fast-oscillating 
perturbations. So we again introduce a slow variable, the angle 8 = n9 — vt. It is also convenient 
to introduce the normalized action I = I/n. The dynamics of the slow variables {1,9} is then 
governed by the auxiliary Hamiltonian H = nH — vl (as may readily be checked by direct 
substitution into the Hamiltonian equations of motion) . Hence 




j = J_H i= 9_H 
89 dl 
9 = n8-vt, I = I/n, 
H = H S + V f , 

H s = H S {I,9) = e ^io 2 J n (kp(I))cos(e), 

V, = V f (I,~9,t) = ^ 2 E ^(M/))cos {^8- (l - ^) vt 

Thus, H s is an autonomous Hamiltonian that determines the main features of the motion of 
{1,9}, while Vf plays the role of a fast-oscillating perturbation. 

It is straightforward to show that the autonomous Hamiltonian H s possesses a single infinite 
grid- like separatrix corresponding to the zero value of H s (Fig. [SJ left panel). The vertical fil- 
aments of the grid correspond to 9 being equal to odd multiples of ir/2 while the horizontal 
filaments correspond to zeros of the relevant Bessel function, 

separatrix : H s = : (20) 
fl = (2j + l)|, j = 0,±l,±2,..., 

I = Ii, i = 0,1,2,..., 
J n (kp{h)) = 0. 

Note that the grid-like separatrix does not depend on the amplitude of the original perturbation. 
Rather its form is an inherent property of the harmonic oscillator driven by the resonant plane 
wave. 

The fast-oscillating term Vf replaces this grid- like separatrix by the narrow chaotic layer. If 
the separatrix (20) is represented in the Poincare section p — q, it takes precisely the cobweb 
form shown schematically in the right-hand panel of Fig. Thus we have achieved the primary 
goal of this subsection, to explain the onset of the cobweb-like stochastic web. 
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Figure 6. Schematic diagram showing the width of the chaotic layer of the web. 




Figure 7. Typical dependence (schematic) of the width of a separatrix chaotic layer on the ratio between the frequency of 
perturbation and the frequency of the small-amplitude oscillation of the unperturbed system. 



2.2.2 Width of the cobweb 

Transport through the web is affected, not only by the shape of the web's skeleton, but also 
by its width, i.e. the width of the chaotic layer (Fig. [6]). An exact calculation of the width is 
a complicated task that we will not undertake here. Rather, we will make a rough estimate 
sufficient to lead us to definite qualitative conclusions. 

Before doing so, we make a general comment about the width in the case of a ID system with a 
separatrix that is being perturbed by a time-periodic perturbation. The width depends strongly 
on the ratio u>f between the frequency of the perturbation w pe rturbation and the frequency of small- 
amplitude eigenoscillations u; unpe rturbed ■ A schematic representation of the typical dependence is 
shown in Fig. [7J This figure will be discussed in more detail in Section O In the present context it 
is sufficient to emphasize that, if iVf is large, then the width of the chaotic layer is exponentially 
narrow. 

Let us now turn to the case of the web. As seen from (19), the characteristic frequency of the 
perturbation Vf is ~ ojq. On the other hand, the unperturbed Hamiltonian H s is proportional 
to the small parameter e. Therefore, even without a careful analysis of its oscillations, it is clear 
that the frequency of oscillation in any cell of its grid-like separatrix is also small. Thus, we 
conclude that the ratio ojf between the perturbation frequency and the eigenfrequency is large, 
so that the width of the layer should be exponentially small. This conclusion is confirmed both 
by careful theoretical analysis and by numerical simulations [H, 27]. 

Moreover, the analysis of oscillations near the elliptic points inside the cells of the separatrix 
of H s shows that, for cells far from the centre, the frequency of oscillation possesses the following 
property [l|, H3] 



^unperturbed (X 



J3/4' 



(21) 



i.e. it decreases as the distance from the centre increases. Conversely, the ratio ut increases. 
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Figure 8. The separatrix for inexact resonance: the single grid-like separatrix of Fig. [5]is replaced by a set of separatrices 
that are distinctly separated from each other. The parameters used to compute the separatrix were: Ao; = —0.001 ^ (see 
Eq. (22)), e = 0.573 and n = k = 1. 

This means that the width of the layer decreases exponentially quickly as the distance from the 
centre of the web increases. This conclusion is confirmed by Fig. 0] above: even for the moderate 
e used in this case, the width of the layer markedly decreases as the distance from the centre 
grows. 

2.2.3 Inexact resonance 

Natural questions to ask in relation to the cobweb are: what happens if the oscillator differs 
slightly from an ideal harmonic oscillator; and what happens if the resonance is inexact? The 
answers were given by Chernikov et al. [23] (see also [l[). They found that the effects of anhar- 
monicity and inexact resonance are in fact similar. So in what follows we shall, for the sake of 
brevity, consider only the inexactness of the resonance: 

v = nu>o + Acj, Alo <C uq. (22) 

In this case, the autonomous resonance Hamiltonian reads as (cf. (19)) 

H s = Aul + ^up n {kp(I)) cos(0) (23) 

As before, there are saddle points corresponding to different values of /, namely different roots 
of the equation J n (kp(I)) = 0. But this means that, unlike the resonance case (Aw = 0), the 
values of H s at the saddles corresponding to different / are themselves different. This means 
that the single grid-like separatrix splits into infinitely many different separatrices (Fig. [8]) . 

In order for at least two lowest separatrices to be connected, allowing chaotic transport within 
a unified structure (a web of a finite size, in the p — q plane), the width of the chaotic layer 
should be more than, or of the order of, the difference in H s between the two lowest separatrices: 

AH > \H S (I 2 ) - H s (h)\ ~ 2vr|Aw|. (24) 

Because AH is exponentially narrow, as discussed in Section [2l2.2 above, the inequality (24) 
means that the stochastic web may be formed only if the perturbation frequency lies in an 
exponentially small vicinity of the resonance. 

2.2.4 Uniform web 

As already demonstrated above, the cobweb cannot in practice provide transport to arbitrarily 
large energies because of the exponentially fast decrease in the width of the chaotic layer with 
distance from the centre of the web. This limitation is overcome in another type of the stochastic 
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web, called the uniform web [30[] (see also [l]). Here, instead of being perturbed by a plane wave, 
the harmonic oscillator is perturbed by short kicks that are periodic in space and time such that 
the kick frequency is equal to the eigenfrequency of the oscillator or to one of its multiples: 

oo 

q + q = — esin(feg) 5{t — nT), (25) 

n=— oo 

T=— , i/ = 1,2,3,... 
v 



The web then covers the whole phase space uniformly, as shown in Fig. [9j 

We note however that the width of the chaotic layer is still exponentially small if the amplitude 
of the perturbation is small [l|, H^j . 



3 Modified stochastic webs 



It is clear from the above discussion that a serious limitation affecting transport through any 
chaotic web is the exponential narrowness of the web's chaotic layer, which leads to exponentially 
slow transport. Soskin et al. [H, 31, 32] recently suggested a way of overcoming this problem by 
making a subtle modification of the webs leading, in turn, to exponential growth in the width of 
the chaotic layer. We shall demonstrate this idea on our example of the cobweb, both because 
it is relevant to the application to the semiconductor SLs and because, in this case, it also leads 
to a dramatic increase in the size of the web. 



3.1 Exact resonance case 



We have found that there is an inherent limitation in the size of the cobweb. It does not relate 
to the inevitably finite time of numerical simulations, which places a practical limit on the 
distance over which the transport can be followed, but is characteristic of the cobweb itself. Our 
numerical simulations show (Fig. [TO]) that, for the given parameters, the inner two-and-a-half 
loops of the web are distinctly separated from the adjacent outer one-and-a-half loops by regular 
trajectories. This might possibly be accounted for theoretically by consideration of higher-order 
approximations of the averaging method [33J. We may speculate that such an approach could 
show that, instead of a single infinite cobweb skeleton, there are many separate separatrices (of 
the one-and-a-half loop shape) lying closely together, but that they might then coalesce due 
to the chaotic layers dressing them as a result of the perturbation. Because the width of the 
layer decreases exponentially fast with increasing distance from the centre, this would mean 
that coalescence would occur only within a few inner loops. Just this is observed in Fig. 1 101 even 
despite that e is moderate rather than small. 
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Figure 10. Left figure. Poincare section for the system q + q = cos(q — t). The saddles (marked by green circles) have 
been found numerically. Four inner saddles belong to one and the same chaotic trajectory (shown in black) which forms 
two-and-a-half inner loops of the stochastic web. Two remaining (outer) saddles generate another (shown in red) chaotic 
trajectory which covers a very thin chaotic layer and is distinctly separated from the black chaotic trajectory. Right figure. 
The area within the dashed magenta rectangle of the left figure is shown on a larger scale. Apart from the black and red 
chaotic trajectories, we show (in magenta and blue) examples of regular trajectories (corresponding to invariant tori) lying 
in between the chaotic trajectories. 



One may reasonably ask: Is there any subtle way to substantially increase the size of the web 
and to enhance transport through it? 

In order to answer this question, let us recall the reason for the exponential narrowness of the 
chaotic layer. It follows from Fig. [7]that it is attributable to the frequency of the perturbation Vf 
being much higher than the eigenfrequency of the unperturbed resonant Hamiltonian H s . It is 
clear from Fig. [7] that the width of the layer would be much larger if we could manage to modify 
the original system in such a way that a new perturbation of the resonance Hamiltonian had a 
component whose frequency was of the order of, or less than, the eigenfrequency of the resonance 
Hamiltonian H s . In fact, this may readily be accomplished in at least two different ways: (i) one 
can add to the original plane wave a small perturbation of the slightly shifted frequency (it can 
itself be e.g. a plane wave); (ii) one can modulate weakly the angle of the original plane wave at 
a low frequency. We demonstrate below only the second option (it will be especially convenient 
for realization of the phenomena in SLs, as shown in Section |4] below). 

We therefore consider the modified system (cf. the original Eq. (15)): 



q + u)^q = e^y- sm(kq — vt — /isin(fii)), 



(26) 



nu , 



n 



1, 2, 3, 



h <C 1, < ^unperturbed ~ ~jz[i 



Of course, the latter inequality cannot be satisfied for an arbitrarily large /, but it can be true 
for a sufficiently high value of I which greatly exceeds the original cobweb size in terms of I. 

Repeating the same procedure used above in the derivation of Eq. (19), i.e. transforming to 
action-angle variables, introducing the slow angle 9 and the auxiliary Hamiltonian H = nH — vl 
which governs the dynamics of {1 — 9}, we can derive: 



H= (modified) 



(27) 



^(modified) ^ ^ + h ™ u lj n ( kp (t)) sin (£) sin (Ot) 
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q/2n q/2n 

Figure 11. Poincare section for a trajectory of the system (28) with initial state q = 0.1, q = (at instants t n = nT where 
T = 2tt/0.02 is the period of the modulation and n = 1,2, 3, ...600000) for ft = (left panel) and h = 0.1 (right panel). A 
symplectic integration scheme of the fourth order is used, with an integration step ti n t = 4 qJ[jq ~ 1-57 X 10 — 4 , so that the 
inaccuracy at each step is of the order of t| nt pa xl0 — 19 . The left panel corresponds to the conventional case considered 
in f!l 1171. \2T\i . The right panel demonstrates that the modulation, although weak, greatly enlarges the web size (note the 
different axes scales), thereby greatly enhancing the chaotic transport. The inset in the top right hand-corner plots the 
left-hand panel on the same scale, thereby illustrating the dramatic extent of this enlargement. 




(in the derivation, we took into account in particular the smallness of h). 

Unlike the original autonomous slow resonance Hamiltonian H s = H S (I,9), the modified slow 
part of the Hamiltonian, i.e. ^ modlfied ) ^ depends on time: it contains a term oc h which oscillates 
at a low frequency f2. It is this slowly oscillating additional term (rather than the former fast- 
oscillating perturbation term Vf) that now determines the width of the layer: the width is 
moderately small (due to the smallness of h and e), rather than exponentially small as in the 
original setup. This exponential growth in the width of the layer gives rise to substantial growth 
in the size of the cobweb. 

To illustrate the above ideas, we use the following example: 

q + q = 0.1 sin[15<?-4t-/isin(0.02i)]. (28) 

For h = 0, this coincides with the conventional cobweb example developed in [H, [l7, 27]. 

Comparison of the left and right panels of Fig. [TTT corresponding to h = and h = 0.1 
respectively, reveals a 6-fold increase in the size of the web in terms of q and p = dq/dt: 



n q! p « 6. 



(29) 
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Figure 13. (a) The potential U(q) = (0.2 — sin(g)) 2 /2, (b) the separatrices in the phase space, and (c) the frequency of 
oscillation as a function of energy uj{E) for the autonomous potential system Ho(p,q, ) = p 2 /2 + U(q). 
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Figure 14. The bifurcation diagram in the plane of amplitude and frequency of perturbation for the system H = Ho + 
hqcos(u>ft) with Ho as in Fig. 1131 The area of {h,uif} for which there is a global chaos between the separatrices is shaded. 
The lower boundary of the shaded area therefore corresponds to the onset of global chaos, representing the function /i cr (u/). 

We emphasize that the modulation giving rise to this substantial increase is actually very small: 
its amplitude of 0.1 is about 60 times smaller than 2tt which is the relevant scale for the angle. 
The corresponding increase of the size in terms of energy is proportional to the square of n qiP : 



n E 



n 



q,p 



36. 



(30) 



Fig- HI shows this explicitly and, in addition, demonstrates that the mode of transport is signif- 
icantly changed. 



3.2 Inexact resonance 

Our idea of an additional small modulation of the angle of the plane wave is equally fruitful 
in the case of an inexact resonance. The frequency band (around the resonance) in which the 
web-like structure is formed may grow exponentially: instead of the exponentially narrow band 
found in the absence of modulation, we may have a moderately narrow band. 

Moreover, there is a nontrivial spectral dependence of this growth: it reflects a universal 
mechanism for facilitation of the onset of chaos between adjacent separatrices, discovered recently 
by Soskin, Mannella and Yevtushenko [34 ]. To explain this mechanism, we use their example: 
it is a potential system with a spatially periodic potential possessing two barriers of different 
height within one period (Fig. 113( a)). Naturally, there are two kinds of separatrices (Fig. 113(b)). 
It was shown j34| that the frequency of oscillation u as a function of energy E possesses a local 
maximum uj m between the separatrices and, moreover, u> is close to oj m over most of the inter- 
separatrix energy range (Fig. [TBT c)). The latter property is valid for any system with two or 
more separatrices and is particularly important in the present context. 

If we perturb the system with a time-periodic perturbation of frequency slightly lower than 
uj m then, due to the flatness of uj(E) over most of the inter-separatrix range of E, two nonlinear 
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Figure 15. Frequency of oscillation in the autonomous Hamiltonian system (23) with parameters as in (31), as a function 
of the corresponding energy E = H s . 
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Figure 16. The spectral dependence (i.e. dependence on f2) of the critical amplitude h of the modulation required for initial 
web formation i.e. for chaotic connection of the first two separatrices of H s (23) to occur. Note the logarithmic vertical 
scale. 

resonances arise that are very wide in terms of energy. Even a rather small amplitude of per- 
turbation may be sufficient for these nonlinear resonances to overlap with each other and with 
the separatrix chaotic layers, thus connecting the latter by the chaotic transport. This has been 
confirmed both theoretically and in numerical simulations. Consequently, a critical perturbation 
amplitude h cr is required for chaotic transport between the separatrix chaotic layers (which may 
be considered as the onset of global chaos between them). As a function of the perturbation 
frequency Uf, it possesses a deep minimu at a frequency approximately equal to u m . This is 
not only seen in the simulations (Fig. [Ti|) but is also well described by the theory [33] • 

The situation is similar for modulation-assisted formation of the web in the case of inexact 
resonance between the plane wave frequency and that of the oscillator. To demonstrate this, we 
use the following example (the parameters correspond to those used in experiments on semicon- 
ductor SLs): 

q + q = esinfg — vt — /isin(Oi)], (31) 
v = 1.02292, e = 0.573. 

For h = 0, a stochastic web is not formed because the Au; = v — 1 k. 0.023 is too large for the 



1 If the perturbation is parametric rather than additive, then the deepest minimum may occur at some multiple of u> m rather 
than at u m itself [34]. 
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chaotic connection of any separatrices of H s (23) to occur. We have calculated numerically the 
two lowest separatrices in the plane 1 — 9 (cf. Fig. [8]), and then obtained the frequency uo of 
oscillation of / (or, equivalently, of the shift by 2ir of 6) as a function of the auxiliary energy 
E = H s : see Fig. [15j There is a local maximum that is clearly similar to that in Fig. [T3Tc). 

Then we switch on the modulation of the wave angle and, for each given Q, increase h gradu- 
ally, until the web is formed i.e. until chaotic connection occurs between the first two separatrices 
H s . This may be considered as the formation of the web. The spectral dependence of the corre- 
sponding critical amplitude is shown in Fig. [T6j Similar to Fig. [bil it exhibits a deep minimum 
(note the logarithmic scale) at a frequency which is a little smaller than the local maximum of 
the dependence u(E): cf. Fig. 16. 



4 Semiconductor superlattices in electric and magnetic fields 

An application of the stochastic cobweb in nanoscience was recently identified and discussed in 
a series of publications by researchers from the University of Nottingham [H, 0, Q, @, [13] • They 
considered quantum electron transport in nanometre-scale ID semiconductor SLs subject to a 
constant electric field along the SL axis and to a constant magnetic field (Fig. [TTla.b)). The 
spatial periodicity of the SL layers gives rise to minibands for the electrons (Fig. I17tc)). In the 
tight-binding approximation, the electron's energy E as a function of its momentum p in the 
lowest miniband is given by [1, 0] 



A[1-cob(M/E)] + gg+rf (32) 



where x is oriented along the SL axis, A is the miniband width, d is the SL period, and m* is 
the electron effective mass for motion in the transverse (i.e. y-z) direction. 

Thus, the quasi-classical motion of an electron of charge e in an electric field F and a magnetic 
field B can be described by: 



f t =-e{F+[V f E(p)xB}}. (33) 



It was shown in [5|] that, for a constant electric field along the SL axis F = (-F, 0, 0) and constant 
magnetic field with a given angle to the axis B = (B cos(#), 0, B sin(#)), the dynamics of the z- 
component of momentum p z reduces to the equation of motion of an auxiliary harmonic oscillator 
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subject to a plane wave i.e. to the equation considered in the previous sections^!!: 

p z + Lolp z = e-r^ sin(/cp 2 -ut + cp), (34) 

LO = UJ C COS(0), U) c =—, 

m* 

eFd 

iy = uj B , w B = — — , 
n 

_ Am* d 2 tan 2 (g) 
e ~ h 2 ' 

dtan(0) 
k = —^> 

<t> = 7T+ -Ap x {t = 0) +p z (t = 0) tan(0)]. 

We emphasize that, despite its classical appearance, Eq. (34) has an inherently quantum origin: 
most of the parameters contain Planck's constant K. 

The dynamics of the system is fully determined by the dynamics of p z . Fig. [18] shows how the 
trajectory of an electron in the x-z plane changes with the angle of the magnetic field. 

At 9 = 0, the plane wave has zero amplitude and the motion along the x- and z-directions 
is separable. Electrons undergo Bloch oscillations along x (due to the presence of the constant 
electric field) and cyclotron motion about B (Fig. flHT a)). The motion is localized. 

Tilting B produces nonlinear coupling of the Bloch and cyclotron motion: as 9 ^ 0, the 
plane wave in (34) acquires a non-zero amplitude. This causes some moderate derealization of 
trajectories (Fig. fTHT b)). The derealization grows very fast (Fig. flST e)) when 9 reaches values 
corresponding to the integer values r = ujb/(u c cos(9)), in other words to the resonance v = uujq. 
This strong derealization in x is a consequence of the onset of the stochastic web for the motion 
of p z (34). 

It is remarkable that the quantum probability density \^(x, z)\ 2 , calculated by solution of the 
Schrodinger equation in the SL model potential subjected to electric and magnetic fields should 
so nicely follow quasi-classical trajectories based on the dynamics of p z (34): see Fig. [I8j 

As shown in [a, 0] , the derealization of the electrons^ strongly affects their drift velocity Vd 
and, as a consequence, the current I and the current- voltage dependence I(V). There are clear 
manifestations of the resonances, both in the theoretical curves Vd(F), I(V), dl/dV (see Fig. 
[TW d). Fig. [T9lf c) and Fig. [T9H ) respectively), and in the experimental curves I(V) and dl/dV 
(Figs. fT9lfa) and (b) respectively). Thus there is clear evidence for stochastic web formation in 
quantum electron transport, providing the basis for a conceptually new method for its control. 
When scattering is included a priori in the semiclassical equations of motion, the stochastic 
web, and stable islands that it enmeshes, evolve into limit cycles. These limit cycles also exhibit 
sharp resonant derealization and their locations in phase space closely reflect the underlying 
web topology [H]. 

Finally, we pose a question: could the modification of the stochastic web discussed in the 



previous section be of use for the SLs? Its seems 2l|, l3ll . |32| ] that this is indeed the case. As 
shown in Section [3] above, modulation of the wave angle results in a large increase of the web 
size. It was noted in [5] that the derealization in x is proportional to the web size in terms of 
the energy of the oscillator in p z , i.e. to E = p z 2 /2 + to 2 p 2 /2 (cf. Fig. [T8|) . The only question is 



1 The only small difference is the presence of a constant shift 4> in the wave angle, but it is inessential. 

2 The motion of electrons in a biased SL with a tilted magnetic field can also be linked to the ultra-fast Fiske effect observed 
for a Josephson junction coupled to an electromagnetic resonator I3. r | . 

1 Seemingly paradoxically, one must take account of scattering in the calculation: if the motion were purely Hamiltonian, 
the position of the electron averaged over time would be constant. 
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Figure 17. (a). Schematic diagram of the SL. Its unit cell comprises two 3.5-nm-thick Ga-As layers (light blue), a 0.3-nm- 
thick InAs layer (dark blue), and a 1-nm-thick AlAs barrier layer (mid blue). The structure contains 14 unit cells, enclosed 
by 50-nm-thick GaAs ohmic contacts (green), (b). Schematic variation of the electronic potential energy with position x 
normal to the layers, for V ^ 0. The quantum wells produce a periodic potential (only two complete unit cells are shown for 
clarity supporting two minibands (blue)). Green areas represent electron gases in the contacts, (c). Energy versus crystal 
momentum dispersion curves for the two minibands. (d). Plots of the drift velocity ly versus F calculated for B = 11T with 
8 = (black curve: arrow marks peak) and 8 = 45° (green curve: arrows mark additional peaks). Reprinted by permission 
from Macmillan Publishers Ltd: [Nature] [H, copyright (2004). 



how the SL should be perturbed in order for the modulation term to appear in the dynamical 
equation for p z 



One suggestion [2l|, [31 



321 ] is that, in a manner similar to the derivation of Eq. 
(34), one can show that the modulation term in the equation for p z appears if an ac-component 
is added to the constant electric field: 



F -> F + F ac cos(uj ac t). 



(35) 



Then, the following modulation term is added in the wave angle in Eq. (34): 



hsm(io ac t), h = — . 

F lOar. 



(36) 



To compare the resulting equation with the example (28) that we studied numerically, we trans- 
form to normalized time 



t — > t = u>ot. 



(37) 
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Figure 18. Electron trajectories and wavefunctions. Grey translucent curves: classical electron trajectories in the x-z plane 
(axes inset) overlaid on corresponding plots of |>I/(:r, z)| 2 (black zero, yellow high) at B = 11T. (a). When 9 = 0, the 
probability distribution is concentrated within the turning points of the classical trajectory, (b). Off resonance, here for 
r = (1 + v5)/4) and 9 = 50°, the trajectories and wave functions extend a little but are still quite localized, (c). On 
resonance, here for r = 1 and 9 = 50°, the wave functions extend across many SL periods, in correspondence with the 
extended classical trajectories. A region of high probability density (yellow peaks) associated with a concentration of orbital 
loops, occurs when the electron is trapped on the first (inner) ring-shaped filament of the web (lower left-hand arrow marks 
the x value corresponding to this ring) and is therefore unable to progress through the SL. But when the electron transfers 
onto the quasi-linear filaments, it shifts rapidly along x, following widely spaced orbital loops (within bracket), which 
correspond to low probability density. The wavefunction is bounded from the right by the second ring-shaped web filament 
(right-hand arrow marks x position corresponding to this ring), which impedes electron flow. Reprinted by permission from 
Macmillan Publishers Ltd: [Nature] copyright (2004). 



Then the equation of motion for p z is: 

d^Pz £ ~ 

-7=5- + Pz = t sin(fcp 2 - vt + 4> + h sin(Qi)), (38) 

dt K 
n _ Woe , _ Fac vJ_UQ 

where all other parameters are as in Eq. (34). 

Thus, if the parameters are similar to those in (28), in particular: h = 0.1, f2 = 0.02, v /ujq = 4, 
then we will have an enlargement in E as found for Eq. (28): tie ~ 36. In order to achieve this, 
we need F ac /F = Ml / (y / ojq) = 1/2000. This means that in order to achieve delocalization of the 
electron by a factor of about 40, we need to add to the constant electric field an ac-component of 
amplitude that is smaller than the constant component by a factor of 2000! We remind the reader 
that the reason for such a dramatic change when an ac-component is added is the exponentially 
strong enhancement of chaotic transport through the stochastic web due to the modulation of 
the wave angle. 

Recently, the effects of stochastic web formation on the high-frequency (GHz-THz) perfor- 



mance of the SLs has been considered 10j, Modulation of the Vd(F) curves, induced by 



stochastic web formation, leads to the formation of multiple propagating electron accumulation 
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Figure 19. Resonant enhancement of current, (a). Experimental I(V) curves measured for B = 11T and 8 = (bottom 
curve) to 90° (top curve) at 5° intervals at a lattice temperature of 4.2K. For 10° < 8 < 55°, each curve contains a region 
of enhanced I beyond V ~ 250mV. (b) Differential conductance plots of the data in (a) reveal strong resonant peaks (red), 
(c). Theoretical I(V) characteristics (for same parameters as in (a)), (d). Differential conductance plots of the traces in (c). 
Curves in (a)-(d) are offset vertically for clarity and those for 8 = 45° are green. In theory and experiment, the resonant 
peaks in G(V) initially shift slightly to higher V as 8 increases, because the enhanced conductance leads to higher electron 
charge density in the SL, which increases F and V . Reprinted by permission from Macmillan Publishers Ltd: [Nature] 
copyright (2004). 

and depletion regions (charge domains), which greatly increase both the strength and frequency 
of the associated temporal current oscillations. Chaos-assisted motion through stochastic webs 
may, therefore, provide a mechanism for controlling the collective dynamics of electrons in SLs 
and, hence, for enhancing their THz performance by using single-particle miniband transport to 
tailor the shape of the Vd(F) curves. 

5 Conclusions 

We have shown that, in general, there is a possibility for energy in a Hamiltonian system to be 
increased from small to rather large values as a result of transport through a stochastic web. 

In a multi-dimensional system, the onset of a stochastic web is a common phenomenon, pre- 
dicted by Arnold in 1964. In the present review, we have been more interested in the low- 
dimensional stochastic webs discovered by Chernikov et al. in the late 1980s. They occur in 
special situations: in a harmonic, or nearly harmonic, oscillator driven by perturbations periodic 
in time and space that are resonant, or nearly resonant, with the oscillator. 

We emphasized that the stochastic cobweb can arise when the oscillator is driven by a weak 
resonant, or nearly resonant, plane wave. The exponentially small width of a strand of the web is 
a characteristic feature of all stochastic webs and it decreases exponentially fast as the distance 
from the centre of the cobweb increases. Moreover there is an inherent limitation on the size of 
the cobweb. Soskin et al. have suggested how to overcome the restriction in size of the cobweb 
and the exponential narrowness of its chaotic layer, just by slightly modifying the system by 
means of a small slow modulation of the angle of the plane wave. 

The model of the stochastic web turned out to be directly relevant to the quantum transport 
of electrons in semiconductor SLs in constant electric and magnetic fields, as demonstrated by 
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Fromhold et al.: the quantum transport dynamics reduces to the model of the harmonic oscillator 
perturbed by a plane wave, where parameters are determined by the values of the electric and 
magnetic fields, by the angle between them, by the period of the SL, by the charge and the 
effective mass of the electron, and by Planck's constant. At certain values of the parameters, 
in particular of the electric field, resonance occurs between the oscillator and the plane wave, 
resulting in the onset of the stochastic cobweb and, consequently, in a strong delocalization of 
the electron which, in turn, increases the current and gives rise to a peak in the dependence of 
the differential conductivity on voltage. 

An addition to the constant electric field of a small slow ac-component results in the slow 
modulation of the plane angle and, therefore, promises to strongly increase the delocalization of 
the electron and to enhance a range of related phenomena. 
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